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Abstract 


In recent years, the occurrence of several destructive floods in the southeastern plain of the 
Caspian Sea, caused by the flood of the Gorganrood River, has caused many casualties and 
damage. Predicting the future morphological conditions of the river and its surroundings is one 
of the essential factors in planning and arranging the coastal plains. Landsat 5, 7, and 8 satellite 
images of 1987, 2002, and 2018, with field studies, software studies, and CA Markov automated 
cells model were used in this research. The accuracy of the modeling was confirmed using random 
ground control points. Moreover, the reaction of the river canal at the flood events was examined 
based on the satellite image of April 3, 2019. The results showed that the most probable changes 
in the study area are in river, plains, and man-made facilities units around Voshmgir Dam Lake. 
Also, the very low hydraulic slope of Gorganrood in the downstream areas of Voshmigir Dam 
prevents the rapid evacuation of floods and creates a reverse flow and flooding of the surrounding 
lands. This phenomenon causes the Gorganrood flood and Qarasoo tributaries to interfere in the 
area between Salagh-Yelghi and Aq Qala. Therefore, increasing the height level of the river 
shores between the tributaries of Gorganrood and Qarasoo can prevent the flow of interference in 
them and reduce the risk of flooding of the Qarasoo around lands. Due to water flooding and flood 
transmission through abandoned canals of Gorganrood in the northern lands of Aq Qala to 
Siminshahr, the construction of an emergency flood drainage canal through abandoned canals can 
be a permanent solution to control floods in these areas. 
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حکیده 


در سال‌های اخیر» وقوع چندین مورد سیل مخرب در دشت جنوب شرقی دریای خزر, ناشی از طغیان آب 
گرگانرود. موجب تلفات و خسارات فراوانی شده است. پیش‌بینی شرایط آینده مورفولوژیک رودخانه و 
محیط‌های پیرامون آن از عوامل ضروری در برنامه‌ریزی و آمایش دشت‌های ساحلی است. در این پژوهش 
از i ait‏ ماهر اتکی a ۷ 8 udo‏ مالا APA e AY‏ هیر اه با وهای دای 
مطالعات نرم‌افزاری و مدل تلفیقی زنجیره مارکوف و سلول‌های خودکار به‌عنوان ipl‏ تحقیق بهره گرفته 
شد. صحت مدل‌سازی با استفاده از نقاط کنترل زمینی تصادفی مورد تأیید قرار گرفت. همچنین نحوه 
عکس العمل QUIS‏ رود در زمان وقوع رخداد سیلاب براساس تصویر ماهواره‌ای سوم آوریل ۲۰۱۹ 
موردبررسی قرار گرفت. نتایج نشان داد بیشترین احتمال تغییرات در منطقه مطالعه در واحدهای رودخانه 
دشت و تأسیسات ساخت انسان پیرامون دریاچه سد وشمگیر وجود دارد. همچنین شیب هیدرولیکی بسیار 
کم گرگانرود در نواحی پایین‌دست سد وشمگین موجب عدم تخلیه سریع سیلاب و ایجاد جریان دبی 
معکوس و طغیان آب به اراضی پیرامون آن اسست. این پدیده موجب تداخل سیلاب گرگانرود و 


سرشاخه‌های قره‌سو در محدوده بین سلاق یلقی و آق‌قلا می‌شود. لذا افزایش تراز ارتفاع کرانه‌های بین 
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سرشاخه‌های گرگانرود و قره‌سو می‌تواند از تداخل جریان در آن‌ها جلوگیری نموده و سبب کاهش 
مخاطره سیل گیری اراضی پیرامون قره‌سو شود. با توجه به طغیان آب و انتقال سیلاب از طریق کانال‌های 
متروک گرگانرود در اراضی شمال آققلا تا سیمین شهر احداث کانال اضطراری تخلیه سیلاب از مسیر 
کانال‌های متروک» می‌تواند یک راهکار دائمی برای کنترل سیل در این مناطق باشد. 

کلیدواژه‌ها: سیلاب. تصاویر ماهواره‌ای» مدل سلول‌های خودکار. مطالعات سنجش ازدور. گرگانرود. 


۱- مقدمه 


مردم در سراسر جهان همه روزه با وقایعی روبرو می‌شوند که منجر به مرگ آسیب. تخریب اموال و اختلال در 
فعالیت‌های روزانه آن‌ها می‌شود. سیل به‌عنوان یکی از مخاطرات طبیعی از دیرباز جوامع روستایی و شهری را مورد 
تهدید و هجوم قرار داده و جوامع بشری فارغ از هرگونه پیشرفت در زمینه‌های مختلف همواره دغدغه‌ی تهدیدات 
سیلاب را به‌عنوان چالشی Soy‏ پیش رو داشته‌اند. استان گلستان در جنوب شرقی دریای خزر» ue‏ سال‌های اخیرء 
شاهد وقوع چندین مورد سیل مخرب بوده که موجب تلفات انسانی و خسارات اقتصادی فراوانی گردیده است. 
رودخانه گرگانرود یکی از مهم‌ترین رودهای سیل‌خیز استان گلستان است. 

فعالیت‌های انسانی با تغییر در رژیم جریان و رسوب. می‌تواند تغییراتی را در پارامترهای مورفولوژی رودخانه به 
وجود آورد (مورایس و همکاران؛ (Y Y‏ بارش‌های شدید و جاری شدن آب از بالادست. ایروبی نشدن 
رودخانه‌ها؛ تغییر در کانال رود و کم شدن فضای آبگذری نسبت به حجم آب ورودی به هنگام بارندگی» از عوامل 
وقوع سیلاب در رودخانه‌های مناطق ساحلی هستند. در این میان اهمیت مطالعات مورفولوژیکی رودخانه‌ها و بررسی 
Se‏ العمل رودخانه در برابر فرایندهای تغییر دهنده برقراری موازنه ضرورتی اجتناب‌ناپذیر اسست (میرزایی‌زاده و 
همکاران ٩‏ ۲۰۱۵). داده‌های ماهواره‌ای یکی از سریع‌ترین و کم‌هزینه‌ترین روش‌های در اختیار محققان» برای تهیه 
نقشه است SL)‏ و مدر ۲۰۰۵). یکی از روش‌های شبیه‌سازی تغییرات در گذر زمان» مدل ترکیبی زنجیره مارکف و 
سلول‌های خودکار* است. با استفاده از این مدل» نسبت تبدیل کاربری‌ها و امکان پیش‌بینی آن‌ها در آینده فراهم 
می‌شود (مولر و VME Po glass‏ بدین ترتیب با استفاده از تغییرات روی داده در گذر زمان می‌توان آینده را پیش‌بینی 
کرد. در سال‌های اخیر رهیافت‌های مفهومی مختلفی از سلول خودکار و تحلیل زنجیره مارکف برای پیش‌بینی 


تغییرات کاربری Os‏ ارائه شده Cal‏ (سوبدی و همکاران ۳ موندال و همکاران" QUY)‏ در تحقیقی تأثیر 
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فیلترهای سازگار با مدل سلول خودکار مارکف را بر نتایج پیش‌بینی تغییر سرزمین بررسی نموده و نتیجه گرفته‌اند که 
ضریب همبستگی O‏ برای فیلترهای ۵6۵ بیشترین تأثیر را در توزیع جغرافیایی/ مکانی داشته است. مدل و همکاران! 
(۲۰۲۰) به بررسی الگوی فرسایش و رسوبگذاری در JUS‏ رودخانه دشت‌های ساحلی پرداخته و نقش بهره‌گیری از 
کانال‌های آنابرانج برای پایدارسازی JUS‏ رودخانه پرداخته‌اند. معصومی و همکاران (۱۶۰۰) به پیش‌بینی تغییرات 
مورفولوژیک رودخانه‌ها در منطقه دشتیاری در کرانه‌های ساحلی دریای عمان پرداخته و عملکرد آن‌ها را طی وقوع 
رخداد سیلاب موردبررسی قرار داده‌اند. ایشان نتیجه گرفته‌اند که سطح مخاطره سیل در پهنه‌های سیلابی منطقه 
ساحلی دشتیاری ory Ub‏ نیازمند عملیات تخصصی کنترل سیلاب از خروجی کوهستان تا dam‏ ساحل است. 
معصومی و همکاران (En)‏ مورفولوژی و فرآیندهای مؤثر در تغیبرات مسیر جریان رودخانه سفیدرود را 
موردبررسی قرار داده و با استفاده از مدل تلفیقی زنجیره مارکوف و clad gle‏ خودکار روند تغییرات مسیر JUS‏ را با 
توجه به عوامل هیدرولیکی رود و تغییرات تراز آب دریای خزر برای سال ۲۰۳۰ پیش‌بینی کرده‌اند. همچنین بر اساس 
نتایج حاصل از مدلسازی» بیان داشتهاند که روند تغیبرات تراز آب دریای خزرب تا سال ۲۰۳۰ به سطح بالاتر از -YV‏ 
متر خواهد رسید. 

با توجه به اهمیت موضوع سیلاب و عملکرد متقابل فرایندهای مؤثر در آن. شناخت این فرایندها اولین گام در 
تحلیل رویداد و اثرات محیطی پیرامون آن است. بر این اساس با توجه به اهمیت موضوع از دیدگاه اقتصادی» 
اجتماعی» زیست‌محبطی و... شناخت تغیبرات رودخانه در طی زمان و فرآیندهای مؤثر از اصئی‌ترین عوامل در 
برنامه‌ریزی‌های مدیریتی برای آینده است. این شناخت با استفاده از تکنیک‌های مدل‌های آماری قابلیت پیش‌بینی در 
بازه‌های زمانی آینده را نیز دارند. با توجه به تعدد وقوع رویداد سیلاب در گرگانرود. شناسایی رفتار رودخانه و 
پیش‌بینی تغییرات مورفولوژی ol‏ می‌تواند کمک شایانی به کنترل سیلاب اراضی پیرامون و برنام‌ریزی در جهت 
عملیاتی نمودن راهکارهایی برای ایمن‌سازی این اراضی در زمان سیلاب نماید. گنجی و همکاران (۱1۰۰) معیارهای 
موثر در سیل‌خیزی رودخانه گرگانرود در محدوده شهرستان آق‌قلا را با استفاده از روش آماری سلسله‌مراتبی AHP‏ 
بررسی نموده‌اند. ایشان نتیجه گرفته‌اند که دیواره‌سازی و تثبیت قوس خارجی پیچان‌رودها در محدوده روستاهای 
shall‏ تا قانقرمه لایروبی رودخانه در نزدیکی de p‏ تغییر کاربری اراضی متجاوز به حریم رودخانه. ایروبی و 
بسترسازی مناسب در محدوده شهری رودخانه گرگانرود و زهکشی و کانال کشی مناسب با شرایط توپوگرافی مناطق 
مرکزی شهرستان از اقداماتی است که به کاهش خطر سیل در شهرستان آق‌قلا کمک می‌نماید. مصفایی و همکاران 
(۱۳۹۹) در نشریه فنی پژوهشکده حفاظت خاک و آبخیزداری. گزارشی از رویداد سیل نوروز ۱۳۹۸ استان گلستان را 


ارائه نموده و عدم مدیریت صحیح مخازن سدهای بالادست و سرعت کم تخلیه جریان در SUS‏ گرگانرود را ازجمله 


1 Medel et al 
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عوامل مؤثر در وقوع این سیل دانسته‌اند. بیت الهی AYIN)‏ در گزارش مطالعاتی» به بررسی دلایل وقوع سیلاب 
۹۸-۷ در استان گلستان پرداخته و نقش عوامل مختلف تأثیرگذار و تأثیرپذیر در این رخداد را مورد تحلیل قرار 
داده است. کاردان و همکاران (Y)‏ این نکته را در مدلسازی دوبعدی سیلاب برای آق‌قلا ذکر کرده‌اند که در سیلاب 
با دوره بازگشت ۱۰۰ ساله» احتمال آب‌گرفتگی و خطر سیلاب بخش عمده‌ای از آق‌قلا را می‌تواند در برگیرد. مزیدی 
و خوش روش (۱۳۹۵) SE‏ تغییر اقلیم بر تعدد رویداد سیل در حوضه گرگانرود را با استفاده از مدل هیدرولوژیک 
بارش- رواناب موردبررسی قرار داده و نتیجه گرفته‌اند که علی‌رغم کاهش در متوسط بارش تعداد بارش‌های حدی 
بیشتر خواهد شد. این امر نشان دهنده افزایش مخاطره سیلاب از نظر تعداد و شدت در آینده خواهد بود. در تحقیق 
حاض با استفاده از داده‌های ماهواره‌ای» میدانی و آماری» مدل تلفیقی زنجیره مارکف و سلول‌های خودکار و نتایج 
تحقیقات گذشته تغییرات مورفولوژی گرگانرود و عوامل محیطی مؤثر بر آن بررسی شده و روند تغییرات و عوامل 
تأثیررگذان برای افق سال ۲۰۳۰ پیش‌بینی شده است. همچنین عملکرد پارامترهای مورفولوژیکی و مؤثر طی رخداد 
سیلاب مورد تحلیل قرار گرفته است. 


۲- محدوده جغرافیایی تحقیق 

این تحقیق بر پیش‌بینی تغییرات مسیر جریان و عملکرد جریان رود در هنگام وقوع سیلاب در رودخانه گرگانرود 
متمرکز شده است. گرگانرود با طولی در حدود ۳۲۵ کیلومتر و حوضه آبریزی به وسعت ۱۰۱۹۷ کیلومترمربع از 
دامنه‌های شمالی البرز شرقی و دامنه‌های غربی ارتفاعات خراسان شمالی سرچشمه می‌گیرد (یوسفی, ۱۳۸۶). جهت 
عمومی جریان آب این رودخانه از شرق به غرب بوده و سرشاخه‌های البرزی آن از جنوب به شمال جریان دارند. 
این رود پس از عبور از شهرهای گنبدک‌اووس و آق‌قلا در غرب خواجه نفس در دلتای گرگانرود به دریای خزر 
می‌ریزد. محدوده مطالعه از ۲۰ کیلومتری بالادست سد وشمگیر تا مصب رودخانه در حدود ۵٩‏ ۳۷۴ تا ۱۵ ۳۷۴ 
شمالی و ۰۰۳ o£?‏ تا 0E? ۵٩‏ شرقی را دربر می‌گیرد (شکل ۱). تاکنون چهار سد در بر روی گرگانرود احداث شده 
که سد وشمگیر و مخازن اصلی و فرعی آن در محدوده این تحقیق واقع شده است. این منطقه بخشی از دشت کم 
شیب شرق دریای خزر را شامل می‌شود. واحدهای چینه شناختی پلیوسن و کواترنری دشت گرگان شامل سازند 
آقچه گیل نهشته‌های مارنی پلیستوسن پیشین» لس‌هاء نهشته‌های گر گانرود. نهشته‌های تبخیری, نهشته‌های برخان‌ها 
و نهشته‌های آبرفتی- دریاچه‌ای می‌باشند. ویژگی‌های رسوب‌شناسی و چینه نگاشتی این واحدها OUS‏ می‌دهد 
فرسایش و نهشته‌گذاری در دشت گرگان حاصل اندرکنش فرایندهای رودخانه‌ای دریایی. دریاچه‌ای و بادی است که 
با تغییرات اقلیمی, فرایندهای زمین ساختی و تغییرات تراز دریای خزر کنترل می‌شده است (قاسمی و همکاران» 
(YA‏ 
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شکل ۱- موقعیت مکانی منطقه مطالعه 


۳- مواد و روش‌ها 

در این پژوهش به‌منظور پیش‌بینی تغییرات مسیر جریان و مورفولوژی رودخانه گرگانرود. از تصاویر ماهواره‌ای؛ 
مطالعات میدانی و مدل‌های نرم‌افزاری استفاده شده است. تصاویر ماهواره‌ای لندست 0 ۷ و ۸ سال‌های ۱۹۸۷۰۲۰۰۳ 
و ۲۰۱۸ با دقت مکانی ۳۰ متر و فصل تصویربرداری مشابه» مبنای مطالعات سنجش زدور می‌باشند. همچنین به‌منظور 
بررسی شیوه عکس العمل SUIS‏ رود در زمان وقوع رخداد سیلاب از تصویر ماهواره‌ای لندست ۸ مربوط به زمان 
سیلاب تاریخ YA‏ اسفند ۱۳۹۷ (تاریخ تصویر سوم آوریل ۲۰۱۹ یا ۱۶ فروردین ۱۳۹۸) استفاده شد (جدول O‏ 
پردازش اولیه تصاویر در نرم‌افزار 53 Envi‏ انجام شد. در این مرحله تصحیحات رادیومتری و اتمسفری به روش 
Quick‏ بر تصاویر ماهواره‌ای اعمال گردید. 

پس از انجام تصحیحات. طبقه‌بندی تصاویر با روش طبقه‌بندی نظارت شده' در نرم‌افزار 10.4.1 ArcGIS‏ 
csi edu ud © io‏ ان و ad‏ کال al steh dle‏ کو cbe‏ کال 
رودخانه» سازه‌های ساخت انسان و سکونتگاه‌هاء اراضی جنگلی» دریا و دشت تفکیک شده‌اند. این نقشه‌ها در هر سه 


دوره تصاویر ماهواره‌ای Angi‏ شده Aul‏ 


1 Maximum Likelihood Classification 


VN‏ جغرافیا و مخاطرات محیطی شمارة سوم 


جدول ۱- مشخصات تصاویر ماهواره‌ای مورداستفاده در پژوهش 


باندهای مورداستفاده تاریخ تصویربرداری مشخصات تصویر 
LANDSAT- LTO05 163-034 1987.06.14 1-2-3-4-5-7‏ 
LANDSAT-LEO07 163-034 2002.07.17 1-2-3-4-5-7‏ 
LANDSAT LCOS8 163-034 2018.07.05 MTL-Multispectral‏ 
LANDSAT LCOS 163-034 2019.04.03 MTL-Multispectral‏ 


به‌منظور پیش‌بینی و مدلسازی نقشه‌های سه دوره به محیط نرم‌افزار Idrisi TerrSet‏ وارد شد. ابتدا براساس 
نقشه‌های سال‌های ۱۹۸۷ و ۲۰۰۲ با استفاده از روش تلفیقی زنجیره‌های مارکوف و مدل سلول‌های خودکار به 
بت سین وو ترات کارمری cs‏ فر ال Y MÀ‏ رداق edd us‏ فار duppi‏ از سیک ها 
احتمالی را بیان می‌کند که از یک وضعیت آغاز شده و به‌صورت پیوسته به وضعیت دیگر تغییر می‌کند (ژانگ و 
همکاران C‏ ۲۰۱۵). مدل تلفیقی زنجیره مارکوف و سلول‌های خودکار عاملی است که تغیبر وضعیت هر سلول را بر 
اساس قانون تعیین وضعیت جدید. مطابق با وضعیت قبلی و وضعیت سلول‌های مجاور نشان می‌دهد. در این مدل 
از فیلتر CA‏ برای توسعه یک فاکتور وزندهی- مجاورت مکانی» برای تغییر وضعیت سلول بر اساس وضعیت 
همسایه‌اش استفاده می‌شود SLL)‏ و همکاران» ۲۰۱۱)؛ بنابراین برای افزايش دقت پیش‌بینی مدل» از این روش 
تلفیقی استفاده شده است. لذا براساس داده‌های خروجی مدل مارکف» نقشه نهایی پیش‌بینی کاربری اراضی با استفاده 
از مدل CA‏ مارکف سال ۲۰۱۸ استخراج گردید و ضریب SUS‏ مربوطه محاسبه شد. ضریب US‏ یک معیار صحت 
سنجی طبقه‌بندی بوده است. سپس خروجی مدل با نقشه حاصل از تصوير ماهواره‌ای سال ۲۰۱۸ برازش داده شد و 
بازه‌های مغایر بین آن‌ها به دست آمد. به‌منظور تعیین دقت طبقه‌بندی و مدلسازی» ۳۷۰ نقطه کنترل زمینی تصادفی 
به‌صورت پیش‌فرض در نرم‌افزار 10.4.1 ArcGIS‏ تعیین شد و با واقعیت مکانی در منطقه موردمطالعه مورد صحت 
سنجی قرار گرفت (جدول ۳) همچنین تراز ارتفاعی نقاطی که منطبق بر کانال رودخانه بودند. جهت ترسیم نیمرخ 
طولی کانال رود برداشت شد. با استفاده از داده‌های ارتفاعی برداشت شده نیمرخ طولی گرگانرود از یلمه سالیان با 
تراز ارتفاعی ۷+ تا مصب رودخانه با تراز ارتفاعی ۲۷- با استفاده از نرم‌افزار Microsoft Excel‏ ترسیم شد. با توجه به 
cds‏ و ey eee‏ مدل با انسفاده از کے هان it ۱۲۰۱۸۵۲۰۸۲ cela Uus‏ پیش کاربری coil‏ و 
مورفولوژی رودخانه برای سال ۲۰۳۰ تهیه شد. این نقشه با برازش خروجی مدل برای سال ۲۰۳۰ و نقشه کاربری 
اراضی سال ۲۰۱۸ به دست آمد. براساس نتایج به دست آمده از این برازش» Geta‏ دارای پتانسیل تغییرات شدید در 


مسیر رودخانه شناسایی شد. به‌منظور بررسی شیوه عکس‌العمل کانال رود در زمان وقوع رخداد سیلاب با استفاده از 


] Zhang et al 
2 Sang et al 
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تصوير ماهواره‌ای لندست ۸ مربوط به زمان سیلاب پهنه‌های متأثر از سیل در نرم‌افزار 10.4.1 ArcGIS‏ در دسا 
تحقیقات بیت الهی (۱۳۹۸» مصفایی و همکاران (۱۳۹۹) و گنجی و همکاران (۱۳۹۹) مورد تدقیق قرار گرفته و تأثیر 
مورفولوژی رودخانه و تغییرات آن در بازه زمانی ۱۹۸۷ تا ۲۰۱۸ و روند آتی آن تا ۲۰۳۰ مورد تحلیل قرار گرفته 


است. 


t‏ نتایج و بحث 

با پردازش و طبقه‌بندی تصاویر در نرم‌افزارهای 104.1 ArcGIS‏ نقشه کاربری اراضی و مورفولوژی رودخانه در 
سه دوره ۱۹۸۷ ۲۰۰۲ و ۲۰۱۸ تهیه شد (Y IC)‏ در این تقشه‌ها واسدهای JUS‏ روشاه SOS Jolt)‏ و بستر 
رود)» سازه‌های ساخت انسان (شامل سکونتگاه‌ها؛ سازه‌ها و ابنیه مختلف)» جنگل» دریا و دشت (پهنه‌های زراعی و 
غیر زراعی) تفکیک شدهاند. مساحت واحدهای طبقه‌بندی شده در هر سه دوره محاسبه شده و مورد مقایسه قرار 


گرفته است (جدول AY‏ 


تصویر ماهواره‌ای 
MANY Yey ۳۰۸‏ 
واحد طبقه‌بندی 
JUS‏ رودخانه ۹Y YYY^VA Y £/A‏ 
تأسیسات ساخت انسان ۱/۷ ۱۳۹۸/۹ enm‏ 
جنگل ۱۳۳۳/۹ ۱۳۹۹/۱ ۱۳۳۳۳۹ 
0V\/0 EAYV MY £AYQYa/£ Cus‏ £42 
Given : Probability of changing to :‏ 


GL. I -D1..2. CL. 3. Gl; 4° Cl. S Cl. 6 


Class 1 : 0.9979 0.0000 0.0021 0.0000 0.0000 0.0000 
Class 2۳: 0.0000 0.8394 0.0505 0.0000 0.1101 0.0000 
Class 3۳: 0.1200 0.0807 0.6746 0.0000 0.1242 0.0005 
Class 4 : 0.0000 0.0000 0.0030 0.9970 0.0000 0.0000 
Class Sam: 0.0003 0.1000 0.1140 0.0000 0.7857 0.0000 
Class 6 : 0.0000 0.0000 0.0180 0.0000 0.0000 0.9820 


شکل ۲- ماتریس خروجی مدل مارکف برای تغییرات سال 1۸ 


vé‏ جغرافیا و مخاطرات محیطی شمارة سوم 


پیش‌بینی و مدل‌سازی» در محیط نرم‌افزار Idrisi TerrSet‏ انجام شد. ابتدا با استفاده از نقشه‌های سال‌های ۱۹۸۷ و 
۲ ماتریس مقادیر احتمال تغییرات کاربری اراضی در سال ۲۰۱۸ بر مبنای زنجیره‌های مارکوف به دست آمد 
(شکل ۳). در این ماتریس مقادیر احتمال تغییرات واحدهای رودخانه (CLS)‏ تأسیسات ساخت انسان (C12)‏ و دشت 
(C13)‏ به یکدیگر مشخص شده است. براین اساس بیشترین احتمال وقوع تغییرات بین واحدهای رودخانه و دشت 
به میزان ۱۳/۵ درصد حاصل شده است (شکل AY‏ 

با توجه به امکان توصیف مکانی در قالب نقشه. از مدل CA.‏ مارکف استفاده شد و نقشه بیش بینی سلول خودکار 
سال ۲۰۱۸ و ضریب GUIS‏ مربوطه استخراج گردید. بر این اساس تشابه مشهودی بین خروجی مدل و نقشه 
طبقه‌بندی تصویر ماهواره‌ای سال ۲۰۱۸ مشاهده شد (شکل (E‏ همچنین ضریب کاپای حاصل شده برای واحدهای 
رودخانه» تأاسیسات ساخت انسان و دشت به ترتیب AAT d. VAR‏ و ۹۸/۲۷ / به دست آمده است. ضریب 
کاپای کلی نیز ۹۸/۷۹ درصد به دست آمده است (جدول (E‏ 

به‌منظور تعیین دقت طبقه‌بندی و مدلسازی» ۲۷۰ نقطه کنترل زمینی تصادفی, با Costly‏ مکانی در منطقه 
موردمطالعه مورد صحت سنجی قرار گرفت. براساس صحت سنجی میدانی. خروجی مدل در ۳۲۱ نقطه صحیح و 
در ٤٤‏ نقطه نادرست بوده است. نتیجه صحت سنجی نشان می‌دهد که در ۸۸/ نقاط کنترل زمینی. خروجی حاصل از 
مدل با واقعیت مطابقت داشته و مطابقت مطلوبی با ضریب کاپای مدل دارد. جزئیات عدم تطابق مشاهدات میدانی در 
جدول Y‏ و ضرایب صحت سنجی مدل در جدول ٤‏ آمده است. براساس نتایج مشخص گردید اصلی‌ترین عدم 
تطابق بین دو واحد دشت و رودخانه به تعداد مجموع alan ٩‏ است. 
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شکل ۳- نقشه طبقه‌بندی و کاربری اراضی و مورفولوژی رودخانه در سه دوره A SAY‏ ۲ و ۲۰۱۸ 


Vo 


M‏ جغرافیا و مخاطرات محیطی شمارة سوم 


شکل -٤‏ نقشه طبقه‌بندی کاربری اراضی و مورفولوژی رودخانه سال ۲۰۱۸ (راست) در مقابل خروجی مدل CA‏ 
مار کف (um)‏ 


با توجه به دقت و صحت خروجی مدل CA‏ مارکف» نقشه پیش‌بینی کاربری اراضی و مورفولوژی رودخانه برای 
سال ۲۰۳۰ با استفاده از نقشه‌های کاربری سال‌های ۲۰۰۲۳ و ۲۰۱۸ تهیه شده است (شکل ۵). سپس سه واحد 
clas eng,‏ مات Oud‏ ور دشست در YA Lus as‏ و ترس ellas s deals, Lp‏ 
محتمل در محیط رودخانه طی بازه زمانی ۱۲ ساله به دست آمد (شکل ۵). 


جدول ۳- جزئیات عدم تطابق مشاهدات در نقاط SAS‏ زمینی با خروجی مدل 


کاربری مشاهداتی کاربری مدلسازی تعداد نقاط 
جنگل دشت ۵ 
دشت رودخانه ۳ 
رودخانه دشت 4 
دشت Lal ese aub‏ ۷ 
دشت E bos‏ 
تاسیسات ساخت انسان دشت ۳ 


با استفاده از این نقشه مناطق دارای پتانسیل تغییرات شدید در مسیر رودخانه شناسایی شد (شکل CV‏ این 
محدوده به طول ۲۳ کیلومتر نواحی پیرامون دریاچه سد وشمگیر از اراضی روستا قزاقلی تا سد را شامل شده و روند 
تغییرات به‌گونه‌ای خواهد بود که واحدهای پیرامون OF‏ دستخوش فرآیندهای فرسایشیی, تغییر پیچانرودی و 
رسوبگذاری خواهند بود. پهنه‌های دارای مخاطره فرسایش و آبگرفتگی (شکل -قرمز رنگ) در جنوب قزاقلی» 
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ساحل جنوبی رود در شمال بی بی شیروان و شمال یارتی قایه مشخص شده‌اند. همچنین ساحل جنوبی رود در 
شمال یاسی e‏ ساحل شمالی رود در شمال بی‌بی شیروان و اراضی حدفاصل بی بی شیروان تا یارتی «B‏ دارای 
قابلیت رسوبگذاری (شکل 7 سبز رنگ) و خروج از آب در پهنه‌های تحت Jb‏ دریاچه سد وشمگیر خواهند بود. 
بر این اساس مشخص می‌شود که بیشترین قابلیت مخاطره فرسایش و آبگرفتگی در اراضی ساحل شمالی دریاچه 
سد وشمگیر وجود داشته و در مقابل» اراضی جنوبی دارای قابلیت رسوبگذاری هستند. مصفایی و همکاران (۱۳۹۹) 
نیز نتیجه گرفته‌اند که حجم رسوبگذاری ذرات گل در محدوده دریاچه سد وشمگیر باتوجه به بافت ریزدانه بار 


رسوبی گرگانرود قابل توجه بوده که این امر موجب تغییرات در موقعیت ساحل دریاچه سد می‌شود. 


جدول ۶- ضرایب صحت سنجی طبقه‌بندی در سه بازه زمانی موردمطالعه 


)/( در سه دوره‎ WIS خطای حذف | دقت کاربر | خطای گماشته | دقت تولیدکننده ضریب‎ 
طبقه‌بندی‎ Joly 
MAY Yey ۸ (A) 9) () C) شده‎ 
رودخانه‎ ۹⁄44 ۸۳/۸ avis ۹۱ 4 AV/AY ۱۳/۷ 
bos aa/£Y Wis) WAE vw ۳۳/۳ WIM YY/N 
جنگل‎ 44/YY aay) ۹۹/۸۳ Ye . ۱.۰ 
دشت‎ ayy AVA 44/YY AY ۱۷ ۸/۷۲۵ \4/Vo 
تأسیسات ساخت انسان‎ ۹۹/۸ ۹5/۸۸ avi‘ ۹۳ ۷ N av 
ضریب صحت کلی‎ AANA ۸۹/۰۹ ANSA 


VA‏ جغرافیا و مخاطرات محیطی شمارة سوم 


| y Uu 

«ind ۲ 
0 A, 2 

b » S. =e. 


شکل ۵- بالا: نقشه طبقه‌بندی کاربری اراضی و مورفولوژی رودخانه eus‏ ۰ خروجی مدل CA‏ مارکف. 
پایین: نقشه پیش‌بینی تغییرات محتمل در محیط رودخانه تا سال YoYo‏ 


54°45'0"E 54°50'0"E 54°55'0"E 
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شکل -٦‏ پیش‌بینی تغیبرات مورفولوژی SUS‏ رودخانه گر گانرود در محدوده تحت تأثیر از دریاچه سد وشمگیر 
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ازجمله آثار سدهای مخزنی روی رودخانه. تجمع و رسوبگذاری آب در مخزن و درنتیجه» کاهش گل آلودگی و 
افزایش ظرفیت نگهداری رسوب جریان در دریاچه سد است. به همین دلیل معمولاً لافاصله در SLL‏ سدهای 
مخزنی» رودخانه دارای ویژگی‌های فرسایندگی بوده و با کف‌کنی بار رسوبی خود را افزایش می‌دهد و پس از رسیدن 
به بخش‌های انتهایی که شیب نیمرخ طولی بسیار اندک می‌باشد. رسوبات ته‌نشین شده و این امر موجب کاهش 
ضریب آبگذری در بخش انتهایی رودخانه می‌نسود (مصفایی و همکاران» ۱۳۹۹). کاهش ضریب آبگذری بخش 
انتهایی رودخانه. عامل مهمی در عدم تخلیه سیلاب از کانال رود به دریای خزر بوده و این پدیده در تمام رویدادهای 
سیلاب موجب آب‌گرفتگی مناطق دشت ساحلی از آق‌قلا تا خواجه نفس شده است. دیگر عامل مؤثر در تغیبر شیب 
هیدرولیکی رودخانه‌های منتهی به دریای خزر» تغییر در تراز آب این درا ات (معصومی و همکاران, ۱۶۰۰). با 
توجه به پیش‌بینی افزایش سسطح تراز آب دریای خزر, طبیعتاً مسطح اساس رودخانه نیز بالاتر خواهد بود. این امر 
موجب کاهش هرچه بیشتر شیب هیدرولیکی رودخانه و تشدید عدم تخلیه آب در هنگام وقوع سیلاب از SUS‏ 
رودخانه گرگانرود خواهد شد. این پدیده موجب سرریز آب از کرانه‌های رودخانه. بخصوص در نواحی نزدیک به 
مصب. در دلتای گرگانرود خواهد شد (شکل ۷ بخش 0( در این زمینه و به‌منظور بررسی نحوه عکس‌العمل کانال 
رود در زمان وقوع رخداد سیلاب. نقشه پهنه‌های تحت نفوذ سیلاب اسفند ۱۳۹۷- فروردین ۱۳۹۸ تهیه شد. سپس 
نقاط حساس در این نقشه با نیمرخ طولی رودخانه در بازه‌ای به طول حدود ۱۳۵ کیلومتر از پایین‌دست سد وشمگیر 
تا مصب» مورد تدقیق قرار گرفت (شکل WV‏ گنجی و همکاران OE)‏ میزان شیب طولی رودخانه گرگانرود در 
محدوده بالادست و پایین‌دست شهر آق‌قلا را ۲ تا Y‏ در هزار تعیین نموده‌اند که با نتایج تحقیق حاضر مطابقت دارد. 
همچنین ایشان مورفولوژی رودخانه در این محدوده را پیچانرود شدید تعیین نموده‌اند. نتایج تحقیق حاضر نشان داده 
متوسط شیب کلی JUS‏ رودخانه گرگانرود در پایین‌دست سد وشمگیر ۰/۰۳ درصد است. متوسط شیب کانال 
رودخانه از خواجه نفس تا سیمین شهر (شکل ۷ بخش ‏ تا ۵) و سلایق یلقی تا انبار آلوم (شکل ۷ بخش ۲ تا ۳ 
برابر با متوسط شیب کلی بازه (۰/۰۳ درصد) و در مقابل شیب متوسط JUS‏ رودخانه در حدفاصل سیمین شهر تا 
GILL‏ یلقی (شکل ۷ بخش ۳ تا (E‏ کمتر از متوسط کلی و به میزان ۰/۰۲ درصد می‌باشد. همچنین متوسط شیب 
کانال گرگانرود از یلمه سالیان تا انبار ell‏ (شکل ۷ بخش ۱ تا ۲) ۰/۰۵ درصد است. نتایج OUS‏ داد که کمترین 
شیب (JUS‏ در حدفاصل سلاق یلقی تا سیمین شهر بوده و این امر موجب کاهش سرعت جریان رو به پایین‌دست 
در این محدوده می‌گردد. لذا در هنگام وقوع سیلاب و جریان با دبی بالاتر از حداکثر توان عبور دهی «JUS‏ در این 
بخش همواره بیشترین احتمال آبگرفتگی وجود دارد. در کنار این مورد کاهش شیب طولی از انبار آلوم به سلاق یلقی 
و ورود سرشاخه فرعی از ارتفاعات البرز به گر گانروده در حدفاصل سلاق بلقی تا آققلا ISS)‏ ۷ فلش آبی SS,‏ 
موجب تشدید آبگرفتگی و ایجاد جریان معکوس از رودخانه به شاخه فرعی می‌شود. تحقیقات مصفایی و همکاران 
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(۱۳۹۹) و بیت الهی (۱۳۹۸) نیز وقوع این پدیده را نشان داده است. در این منطقه با فاصله‌ای حدود ۲۵۰ متر» یکی 
از سرشاخه‌های فرعی قره‌سو به سمت شمال و غرب جریان دارد (شکل ۷ فلش زرد رنگ). وقوع جربان برگشتی 
در گرگانرود موجب سرریز شدن سیلاب شده و جریان سیلاب گرگانرود با فره‌سو پیوسته شده و در مسیر این رود 
موجب سیل گیری اراضی شمال Ade OE S‏ و جنوب بندر ترکمن می‌گردد. چنین پدیده‌ای در سیل اسفند -AV‏ 
فروردین ٩۸‏ سبب تشدید سیلاب در اراضی پیرامون قره‌سو با تراز ارتفاعی حدود ۱۳- متر شده است (شکل (V‏ در 
صورت افزایش تراز اب دریای خزر پس‌زدگی و ایجاد جریان برگشتی» به سطوح ارتفاعی بالاتر نیز کشیده خواهد 
i sa cid‏ هیال ای eas ose sa‏ قفا تساه ناس یات مار ANG deba‏ 
تراز ارتفاع کرانه‌های این دو سرشاخه می‌تواند از تداخل جریان آن‌ها جلوگیری نموده و نقش قابل‌توجهی در کاهش 
مخاطره سیل گیری اراضی پیرامون قره‌سو داشته باشد. بدیهی است این امر نیازمند انجام پژوهش‌های تخصصی در 
زمینه مهندسی رودخانه و سازه‌های حفاظتی مرتبط می‌باشد. نتایج پهنه‌بندی سیلاب فروردین ۱۳۹۸ نشان داد که در 
منطقه سلاق یلقی- آق‌قلاء (شکل V‏ بخش ۳) علاوه بر حروج سیلاب گرگانرود به قره‌سو» بخشی از جریان آب با 
سرریز از کرانه‌های رود. از شمال قره‌سو به سمت گمیش تبه جاری شده است. براساس شواهد تصاویر ماهواره‌ای و 
اقات شلات این عمط بر کال ھا وین sad led‏ کاود انیت کال Vu‏ بیان id‏ 
متروک در این مسیر اراضی گود افتاده‌ای را به وجود آورده‌اند که سیلاب به سمت آن‌ها جریان یافته و از این طریق 


در دشت شمالی گر گانرود پخش می‌شود. 
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و مسیر شاخه فرعی گرگانرود (فلش آبی) و سرشاخه قره‌سو (فلش زرد) 


فعالیت‌های انسانی از قبیل راه‌آهن و راه‌ها ب‌صورت متقاطع با مسیر جریان ذکر شده موجب گسترش سیلاب در 
دشت شمالی گرگانرود می‌گردد. نفوذ سیلاب از بخش شمالی شهرهای WET‏ (شکل ۷ بخش ۳) و سیمین شهر 
(شکل ۷ بخش E‏ در سیل فروردین ۱۳۹۸ در اثر جریان سیلاب از مسیر یاد شده و تحت تأثیر سازه‌های تقاطعی 
صورت گرفته است. جریان سیلاب در این مسیر پس از عبور از اراضی شمال سیمین شهر. به سمت گمیش تبه 
جریان یافته و در شرق گمیش تپه به دو بخش تقسیم شده است. بخشی از جریان از طریق آبراهه‌های زهکشی. از 
شهر گمیش ته به تالاب گمیشان را یافته و بخش دیگر در اراضی پست شرق و شمال شرق این شهر گسترده شده 
است. بیت الهی YA)‏ بیان داشته که شکافتن مسیر راه‌آهن در چندین نقطه در محدوده آق‌قلا در کاهش ارتفاع آب 
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جمع شده در شهر موثر بوده است. نتایج تحقیق حاضر نشان می‌دهد که در سیل فروردین ۱۳۹۸ بخش قابل‌توجهی 
از > Ob‏ سیلاب گرگانرود. از طریق این مسیر تخلیه شده است. لذا لایروبی. مسیر سازی و اصلاح تقاطع راه‌ها و 
راه‌آهن با این مسیر سیلابی و احداث SUIS‏ اضطراری تخلیه سیلاب به اراضی شمالی گمیش تیه می‌تواند یک راهکار 
دائمی برای JES‏ سیل در مناطق پایین‌دست سد وشمگیر باشد. 


0— جمع‌بندی 

وقوع سیل در دشت گرگانرود و تکرار هرساله آن. یکی از مهم‌ترین مخاطرات محیطی در این منطقه است. در 
این میان عدم تخلیه جریان سیلاب از رودخانه تحت تأثیر تغیبرات مورفولوژی SUS‏ و تراز آب دریای خزر dele‏ 
مهم در سیل‌خیزی هستند. با توجه به تغییرات دائمی در تراز آب دریای خزر فرآیندهای رودخانه‌ای و فعالیت‌های 
انسانی» شناسایی تغیبرات رودخانه و پیش‌بینی آن بسیاری ضروری است. تحقیق حاضر در راستای این اهداف 
صورت گرفت و مباحث pd‏ از Ol‏ نتبجه گیری شد. 

- مدل تلفیقی CA‏ مارکوف قابلیت مناسبی از نظر دقت و صحت در پیش‌بینی تغییرات مورفولوژی رودخانه 
گرگانرود نشان داد. 

- بیشترین احتمال تغییرات در منطقه موردمطالعه در واحدهای رودخانه» دشت و تأسیسات ساخت انسان پیرامون 
دریاچه سد وشمگیر از اراضی روستا قزاقلی تا سد وجود دارد. روند این تغییرات در این نواحی cla Sa‏ خواهد 
بود که اراضی جنوبی پیرامون رودخانه دستخوش فرآیندهای فرسایشی. تغییر پیچانرودی و سیل‌گیری خواهند شد. 
در مقابل اراضی شمالی غالا تحت روند رسوبگذاری قرار خواهند داشت. 

-با توجه به شیب هیدرولیکی بسیار کم گرگانرود در نواحی پایین‌دست سد وشمگیر, عدم تخلیه سیلاب و ایجاد 
جریان دبی معکوس, یکی از پدیده‌های مهم در UL‏ آمدن تراز رودخانه و طغیان آب به اراضی پیرامون آن می‌باشد. در 
سیل فروردین ۱۳۹۸ عمق نفوذ این پدیده تا محدوده بین سلاق بلقی و آققلا pli‏ حدود ۱۳- متر) بوده و موجب 
[ols‏ سیلاب گر گانرود و سرشاخه‌های قره‌سو شده است. لذا در صورت افزایش تراز آب دریای خزر این مخاطره 
به اراضی با تراز ارتفاعی SUL‏ نیز اثرگذار خواهد بود. 

- براساس نتایج پهنه‌بندی سیل فروردین ۱۳۹۸ء جریان آب علاوه بر JUS‏ گرگانرود از طریق قره‌سو در جنوب و 
مسیر کانال‌های متروک گرگانرود در شمال آن تخلیه طبیعی شده اسست. انجام پژوهش‌های تخصصی در زمینه 
مهندسی رودخانه و سازه‌های حفاظتی مرتبط در قالب دیوار ساحلی و افزایش تراز ارتفاع کرانه‌های بین سرشاخه‌های 
گرگانرود و قره‌سو می‌تواند از تداخل جریان در آن‌ها جلوگیری نموده و نقش قابل‌توجهی در کاهش مخاطره سیل 


گیری اراضی پیرامون قره‌سو داشته باشد. همچنین احداث کانال اضطراری تخلیه سیلاب از مسیر کانال‌های متروک 
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sbi‏ باشد. 


سپاسگزاری 

بدین‌وسیله از همکاران گروه مهندسی رودخانه و سواحل پژوهشکده حفاظت خاک و آبخیزداری در اجرای 
تحقیق قدردانی می‌گردد. این پژوهش برگرفته از نتایج طرح تحقیقاتی پیش‌بینی روند let‏ مورفولوژی 
رودخانه‌های مناطق ساحلی شمال کشور با کد ۰-۲۹-۲۹-۰۲۳-۹۷۰۵۳۵ پژوهشکده حفاظت خاک و آبخیزداری 
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